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Abstract 

We show that in the transitions ip' — > h c ir° and r\ c — > XcOit the contributions from charmed meson loops are highly suppressed, 
in contrast to various other charmonium decays. We calculate the width of the ip' — > h c ir°, which agrees with the recent BES-III 
data, and predict the width of the n' c — > XcOTr , ^{n' c — ► XcOi" ) = 1.5 ± 0.4 keV. A confirmation of this prediction would also 
provide additional support for a recent analysis of ip' — ► J/ipn°(n), where loops are claimed to play a prominent role. 

PACS numbers: 13.25.Gv, 14.40.Pq, 12.39.Fe 
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I. INTRODUCTION 

Recently, in a variety of calculations it was shown that 
charmed meson loops contribute importantly to the de- 
cays of charmonia with and without isospin breaking (for 
a overview, see [1]). For instance, using an effective La- 
grangian approach (EL A), the intermediate meson loop 
contributions are found to be essential for understand- 
ing the puzzling ip (3770) non-DD decays [2j, |3j. They 
are also important in the J/ip decays into a vector and 
a pseudoscalar mesons Q and in the Ml radiative tran- 
sitions between two charmonia Q. For isospin breaking 
decays, besides ir° — 7) mixing and electromagnetic (e.m.) 
effects, one also expects the mass difference between the 
neutral and charged charmed mesons in the intermediate 
states (i.e. in the meson loops) to play a role. This effect, 
known to be of particular importance near the continuum 
thresholds, was studied in the decays 4> — > W7r ° [E 0h 
J/ip -> (j)r]ir 0,0, and £>* (2317) ->■ D s n° [IM1- 

In Ref. 13J a non-relativistic effective field theory 
(NREFT) was described, that allows one to study the role 
of charmed meson loops in charmonium decays in a sys- 
tematic way. Applied to the reaction ip' — > J/ipw°(r]), the 
formalism revealed that the contribution from charmed 
meson loops is enhanced by a factor of 1/v, with v ~ 0.5 
being the charmed meson velocity [4lJ, compared to the 
tree-level one. In this letter we apply NREFT to the de- 
cays ip' — > h c n° and rj c — > Xcqk ■ We will show that in 
these two decays, the loop contributions are highly sup- 
pressed, and hence the tree-level terms dominate the de- 
cay amplitudes. Testing experimentally the predictions 
that emerge here, especially for the partial decay width 
for i]' c — > XcO 71 " ! would provide a non-trivial test of the 
NREFT and is thus of high importance towards an un- 
derstanding of the properties of charmonia. 

In this context it is instructive to compare the ip' — > 
J I tpir decays to ip' — > h c ir° on the basis of power count- 
ing. Here we follow the reasoning of Ref. [lj|. For the 
former decay, which happens in a p-wave, the tree level 



amplitude scales as m q q, where q denotes the momen- 
tum of the final particles in the ip' rest-frame and a 
quark mass factor appears since the reaction is isospin 
violating. For the same reaction the loops on the other 
hand scale as v 3 /v 4 (m q /v 2 )[qv 2 ] = qm q /v, where the 
the factor v 3 comes from the the non-relativistic inte- 
gral measure and the 1 jv 4 from the two non-relativistic 
two-meson propagators. Further, the term in the round 
brackets emerges, since pulling out a factor of m q , which 
is an energy scale, has to be balanced by a factor which 
characterizes the intrinsic energy, v 2 . In addition, the 
[qv 2 ] term contains the vertex factors from the external 
pion coupling, q, and from the two p-wave vertices in the 
loop. Thus, heavy meson loops are enhanced by a factor 
1/v ~ 2. For the s-wave decay ip' — > h c Tr° on the other 



hand, the tree level scales as 



'Q I 



while the loop here 



scales as v 3 /v (m q /v 2 )[q/MD]' 2 . All factors that appear 
are analogous to those discussed above up to the factor 
q 2 /All). The origin of this is, on the one hand, that the 
pion is produced in a p- wave and, on the other hand, that 
there is now only one p-wave vertex in the loop (the ip' 
decay) giving rise to a momentum factor at that vertex. 
To obtain a non-vanishing result, however, that momen- 
tum has to be proportional to q (c.f. Eq. (fT5")) below). 
The factor l/Mf, is then introduced to match dimen- 
sions. Thus, for the reaction ip' — > h c Tr° loops appear 
to be suppressed kinematically, on the amplitude level, 
by a factor q 2 /(M^v 3 ) ~ 1/30. In case of the reaction 
Vc XcO 7 !" j we find from the same analysis a suppres- 
sion of the loops by factor 1/10 — here the suppression 
is weaker due to a larger phase space. Thus we find on 
the basis of the same power counting that heavy meson 
loops are enhanced compared to the tree level amplitudes 
in the isospin violating p-wave decays like ip' — > J/ipir°, 
while they are strongly suppressed in s-wave decays like 
ip' —> h c TT°. The latter observation allows us to predict 
the partial decay width of rj' c — > XcO 7 ? ■ A confirmation of 
this prediction would at the same time provide a strong 
support for the analysis of Ref. fl3[ . 
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The hcQ-Pi) has been the last charmonium state be- 
low the DD threshold that was confirmed experimen- 
tally; for a comprehensive review of the charmonium 
physics, see Ref. [13]. It was first established in pp an- 
nihilation by the E760 Collaboration at Fermi Lab in 
1992 [lBl - With J PC = this state cannot be pro- 

duced in e + e~ annihilation directly. Furthermore, due 
to the phase space restriction, it cannot be accessed by 
ip' decays into h c rj. Instead, the only open strong de- 
cay to h c is via the isospin-violating ip' — > h c ir° pro- 
cess. As a consequence, this branching ratio is strongly 
suppressed. Recently, the CLEO-c Collaboration with 
24.5 x 10 6 ip' events succeeded in measuring precisely the 
product of two branching ratios B{ip' —> K°h c ) x B(h c — > 
■yrjc) = (4.19 ±0.32 ±0.45) x 10~ 4 in [1^. This combined 
branching ratio was confirmed by the BES-III Collabora- 
tion [ljj with 110 million ip' events in the same channel, 
B(ip' -> n°h c )xB(h c ->• jtj c ) = (4.58±0.40±0.50) x 10~ 4 
Furthermore, they reported the first measurement of the 
absolute value of the branching ratio for the ij/ — > h c Tr° 
as B(ip' -> Tr°h c ) = (8.4 ± 1.3 ± 1.0) x 10" 4 [l8[. Using 
the PDG value for the total width of the ip', T(ip') = 
309 ± 9 keV [H[, the partial width of the ip' h c Tr° is 

Y{ip' h c n°) = 0.26 ± 0.05 keV. (1) 

This experimental progress makes it possible to study 
physics of the h c to some extent of precision. In partic- 
ular, its production in the decay ip' — » h c ir° appears to 
be an ideal channel for investigating the isospin-violating 
mechanisms and the pertinent non-perturbative QCD dy- 
namics. The reaction ip' — s- h c ir° was first studied the- 
oretically more than thirty years ago using k — •// mix- 
ing [201 . In order to get an estimate of the width, the 
authors estimated the coupling of the ip' to the h c and rj 
by assuming it to be equal to the ip' J/ipr] coupling, which 
was extracted from the measured ip' — ¥ J/ipr/ width. As 
a result, they obtained 5 ... 30 keV for the partial decay 
width of the ip' — > h c ir°, which overshoots the measure- 
ment significantly. 

The QCD multipole expansion (QCDME) was applied 
to this problem [2l| - |24j ] , and the following typical expres- 
sion was obtained 

T[ip' -> h c n°] = 0.12— keV (2) 

aE 

where a e and a m are the coupling constants for the color 
electric dipole and magnetic dipole gluon radiation, re- 
spectively. The phenomenological determination of the 
ratio oim/oie has large uncertainties. By taking the ratio 
in a range of 1 ... 3, the partial width is about 0.12 .. . 0.36 
keV, consistent with the experimental result. In contrast, 
a later calculation [25] gives a larger value of 0.84 keV, 
and the estimate by Voloshin gives a much smaller value 
of 15 eV [26] ]. Suffering from a poor knowledge of the 
coupling constants fillip, the QCDME results should 
rather be regarded as an order-of-magnitude estimate. 

In this work, we shall further investigate the isospin 
violation mechanisms of tp' — > h c ir° and its analogue, 



the reaction r]' c —> XcQ^ ■ In Section [HI we will give 
the tree level decay amplitudes by constructing the ef- 
fective chiral Lagrangian. In Section IIIII we present the 
NREFT [l3[ analysis for the heavy meson loops. As 
we will demonstrate the explicit calculation supports the 
scale arguments presented above, that heavy meson loops 
are highly suppressed for the reactions under consider- 
ation. We also checked that the results are consistent 
with a calculation using the ELA d S — details will 
be presented elsewhere 28] . In Section HVl the results 
for the decay widths are given. In particular, the width 
of the r/' c — > Xco 71 " is predicted. Some discussions and a 
summary are given in the last section. 

II. TREE LEVEL CONTRIBUTION 

Since the mass difference between the initial and final 
charmonia is small, the emitted pion is soft. Hence one 
can construct an effective chiral Lagrangian considering 
the charmonia as matter fields. Since the charmonia are 
isoscalar and the pion has isospin one, the transitions vio- 
late isospin symmetry. Isospin breaking has two sources. 
One is the mass difference between the up and down 
quarks, and the other one is of e.m. origin. For the tran- 
sitions considered here, the e.m. effect can be neglected 
(for details, see HH). Defining \ — 2£?o • diag (m u , m^), 
the quark mass difference is contained in the operator \- 
which contains an odd number of pion fields, 

X- = u^xu 1 ~ ux^u, (3) 

where B = \ (0\qq\0) \/F%, is the pion decay con- 
stant in the chiral limit, and u parameterizes the pion 
fields as the Goldstone bosons of the spontaneously bro- 
ken SU(2) L xSU(2) fl , 

' 4)- w 

Using the two-component notation of Ref. [29j , the field 
for the S'-wave charmonia ip' and rj' c reads 

J' = iP'-a + r 1 ' c , (5) 

with ip' and rj' c annihilating the ip' and rj' c states, and a 
the Pauli matrices. The field for the P-wave charmonia 
is 

X 1 = o> (~X% - ^ ijk X k d + ^Xco) + K, (6) 

where x% Xli, XcO and h c annihilate the Xc2, Xa, Xco 
and h c states, respectively. 

The leading order (LO) chiral Lagrangian for the tran- 
sitions ip' — > h c ir° and r( c — > XcO 71 " reads 

£ = \C [<X f ■ 3 J') + (J 'o ■ x f >] (X-) aa + h-c, (7) 
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where C is an unknown coupling constant, (■ • •) is the 
trace in spinor space, the subscript a — u, d is a fla- 
vor index, and repeating indices aa means evaluating the 
trace in flavor space. The spin symmetry is violated due 
to the presence of the Pauli matrices between the two 
heavy quarkonium fields. Note that the Lagrangian was 
first proposed in Ref. (30l | in four-component notation, 
and the coupling constant C here is —4/Bq times the 
one defined in that paper. Working out the traces in 
both the spinor and flavor space, one finds that only the 
two transitions considered in this paper are allowed as 
dictated by conservation of angular momentum, parity 
and charge conjugation invariance, 



C = 



iC (V • hj + V3r,' cX t ) ( X -) aa + h.c, (8) 



and, after taking into account the ir° — rj mixing which 
contributes to the isospin breaking transitions, 



(X-)aa = 6i^(m d - m u )n a + 

-Fir 



(9) 



where tt° = ir a + e 7T ° v r ] is the physical pion field with e^o^ 
being the tt° — r\ mixing angle. In the above equation, we 
have neglected the multi-pion terms. Now it is easy to 
write out the tree level amplitudes, 

G 



M(tf -> Kit ) = y/M^Mh—CW) ■ e(h c ) 

xB (m d - m u ), 
M(r)' c -> Xcott ) = ^jM v ,M Xc0 ^-CB (m d - m u %10) 

Note that the above amplitudes were multiplied by a fac- 
tor y/MfMi, with Mjffj being the mass of the initial 
(final) charmonium, to account for the non-relativistic 
normalization of the heavy fields in the Lagrangian. 

III. CHARMED MESON LOOPS 

In this section we investigate the contribution of 
charmed meson loops to the two decays. The relevant 
diagrams for the reaction tp' — > h c ir° with neutral in- 
termediate states are shown in Fig. [TJ The coupling of 
pion to the charmed mesons is described by heavy me- 
son chiral perturbation theory p3iT[33j (for a review, see 
Ref. |34[). The fields for the pseudoscalar and vector 
charmed mesons in the same spin multiplet can be writ- 
ten as H a = V a -<7 + P a with V a and P a denoting the vec- 
tor and pseudoscalar charmed mesons [2{|, respectively, 
where a is the flavor index with {P u ,Pd} = {D°,D + } 



and similar for the vector mesons. The LO chiral effec- 
tive Lagrangian for the axial coupling is [29| 



(ID 



where the axial current is u = 
g the pertinent coupling constant. 




FIG. 1: Diagrams of neutral intermediate meson loops con- 
tributing to ip' — > h c ir°. Their charge conjugate diagrams are 
implied. 



The LO Lagrangian for the coupling of the S- or 
P-wave charmonium fields to the charmed and anti- 
charmed mesons can be constructed respecting parity, 
charge conjugation and spin symmetry. The one for the 
if)' and T]' c reads [l3| 



(12) 



where A dB = A{dB) - (dA)B, and //,. -V a ■ a + P a 
is the field for anti-charmed mesons [35|. The LO La- 
grangian for the P-wave charmonia spin multiplet is (35j 



(13) 



These Lagrangians were introduced in Ref. [36| in four- 
component notation. The values of the coupling con- 
stants gi and 172 in that paper are half of those introduced 
here. 

Using these Lagrangians, one can work out the decay 
amplitudes. For simplicity, we focus on diagram (a) in 
Fig. [T] — the analysis of the other diagrams is analogous. 
The amplitude for ip' — > h c ir° from the diagram (a) reads 
(for the full amplitudes and calculation details, we refer 
to Ref. [H) 



J 



M(i/>' -)• fc c 7r°) w 



fir 



e(h c )q- eW)R \l^(q, D\ D°, D*°) - D ± ,D ± ,D* ± ) 



(14) 
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where q is the pion three-momentum in the rest-frame of 
the initial charmonium, R = w 'M^/M^Mp M o * is in- 
cluded to account for the non-relativistic normalization 
of the heavy fields, and we have labeled the loop func- 
tion I^'(q) by the intermediate charmed mesons. It is 



obvious that the isospin breaking comes from the neutral 
and charged meson mass differences. The loop function, 
evaluated in the rest-frame of the decaying particle, is 
defined by 



8mi 



dH 



8m im2 m 3 J (2tt) 4 (jo_ ri (|jM fpo_jo_ ra (|j|)) 6o_ g o_ T3 (|f_^) 

i r d 3 i r 



m 2 m 3 J (2tt)3 ^_ ra (|i|)_ Tl (|j|)) (E f -T a (\l\)-T 3 (\T-%) 



(15) 



where Tj(p) = p 2 /2mi denotes the kinetic energy for 
the charmed mesons with masses mi,m 2 and m 3 , Ei — 
Mi — mi — ?7i2 and £y = Mf — rri2 — m 3 — E n denote the 
energies available for the first (before pion emission) and 
second (after pion emission) two heavy meson intermedi- 
ate state. For diagram (a), mi = mi — Md, rn 3 — Md*, 
Mi = Mm, and M f = M hc . The loop function is con- 
vergent [H, [H|. Defining c = 2/Ji 2 fri2 and c' — 2/J23623, 



f/(r 



with Hij 

bi2 = m i + m 2 - Mm and 623 = m 2 + m 3 
the loop function can be approximated by 



being the reduced mass, 

E n — Mm, 



I (1 \q)=N 



M23 



2^~c 



m 3 



(16) 



with N = /xi2A*23/(167rmim2TO3), where terms of order 
q 2 /c' and higher have been neglected. The approxima- 
tion is reasonable because for either of the two decays 
considered here, the pion momentum is small and fulfills 
q 2 < c'. 

We can now have another, more refined look, at the or- 
der of magnitude estimate for the loop function. Since yfc 
and \fd are approximately the momenta of the charmed 
mesons in the loop, we count them as Mj^v with v be- 
ing the velocity of the charmed mesons. It follows that 
&12 ~ &23 ~ Mdv 2 . For an order-of-magnitude estimate, 
one may neglect the difference between c and d , and 
denoting them by 2[ib. Then, I (1) (<?) ~ N/(4^/2fIb). De- 
noting the mass difference between the charged and neu- 
tral charmed mesons by 5, we have fi c = fi n + 5/2 and 
b c = b n + 25, where the lower-index n(c) means neutral 
(charged). The amplitude for diagram (a) scales as 



(a)" 



N x q 



5v 



b' 2 



?! 2 



AT N X 1 * 

,N ^n 5 ^Mj' 



(17) 



where N {a} = 2(g/F n ) 9l g 2 Rq ■ e{h c )q ■ e(^')/q 2 . Thus, 
we confirm the parametric behavior derived in the intro- 
duction on the basis of the NREFT. For the transition 
from the ip' to the h Cl the pion momentum in the ip' rest- 
frame is q = qi = 86 MeV, and hence q 2 /M|, ~ 2 x 10~ 3 . 
Taking into account that the velocity v may be roughly 
estimated as y/[2Mf, - (M^> + M h P)/2\/M b ~ 0.4 with 
Mjj being the averaged charmed meson mass, the dimen- 
sionless factor 



v 3 M 



0.03 



(18) 



D 



produces a significant suppression compared to the tree 
level amplitude. A similar though more moderate sup- 
pression happens in case of the rj' c — > XcO 71 " also. For this 



decay, the momentum of the pion is q = q% 
and the suppression factor is 



1 



■J 3 M 2 D 



0.1. 



171 MeV, 



(19) 



The numerical results support the above power count- 
ing argument. If we only consider the contribution from 
the charmed-meson loops, the widths of the ip' — > h c Tr° 
and rj' c — > XcO^° are 



T(ip' -> h c ir°)i oop - 

nv'c -> XcO^°)ioo P = 1.0 x 10^ g 2 g' 2 2 keV, (20) 



= 2.1 x 10- 7 ff 2 ff2 2 keV, 

5„2 / 2 



where the 7r° — r\ mixing has been taken into account, and 
the values of g\ and g' 2 are given in units of GeV -1 / 2 and 
GeV~ 3 / 2 , respectively. We checked that the EL A gives 
similar results, which confirms our analysis. 

One may ask if the values of the coupling constants 
are so large that the suppression gets invalidated. In 
fact, because all the charmonia considered here are be- 
low the DD threshold, the couplings cannot be extracted 
directly from the decay widths. However, one may get 
a feeling about their values from other sources or from 
model calculations. Assuming the coupling of the ip' to 
the charmed mesons has similar strength as the one of the 
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J/ip, we obtain g' 2 ~ 2 GeV~ 3 / 2 from the ^' — ► J/ipir°(ri) 
where the charmed meson loops dominate 13]. Values of 
the same order of magnitude were obtained from various 
model calculations, see e.g. Refs. [3w3o |. In Ref. [36| . 
the authors estimated g\ using vector meson dominance, 
which gives g± = —4.2 GeV -1 / 2 . Using these values, 
the resulting width F(ip' -> h c ir°)i oop ~lx 10~ 5 keV. 
It is smaller by 4 orders of magnitude compared to the 
BES-III measurement, and it confirms our power count- 
ing estimate presented before. 

In addition, even without any assumption on the cou- 
pling constants, from our analysis we can predict 



^(Vc XcQ7r°)loop 

r(^ -> h c n% oop 



= 48 



(21) 



As we will see, this ratio as derived solely from the loop 
contributions, is much larger than the corresponding one 
derived from including the tree level amplitudes only. 
Note, however, that because v ~ 0.5, there might be 
sizeable corrections to this result. Thus, for the men- 
tioned decays it can be tested experimentally, if there is 
a dominance from the loops or from the tree level contri- 
bution. 



the total width of the rf c , r(r)' c ) = 14 ± 7 MeV [1J, the 
branching fraction of the isospin-breaking transition is 



B (»7c "> Xcott ) = (1.1 ±0.6) x 10" 



(25) 



Note, as a contrast, if the reactions were dominated 
by heavy meson loops, the predicted branching fraction 
would be larger by a factor of about 5. The predic- 
tion, and therefore the dynamics underlying the decays, 
is testable with PANDA at FAIR [Hj]. 

Furthermore, one may check if the tree-level contri- 
bution gives the right order of magnitude of the decay 
width T (ip' — > h c Tr°) through dimensional analysis. The 
tree-level amplitudes of Eq. (fTU)) are proportional to the 



dimensionless factor ^JMfMiC . Because the spin sym- 
metry is violated as can be seen from the presence of the 
Pauli matrices in Eq. ([?])■ one may write 



rn € 



(26) 



with the dimensionless parameter C being a number of 
natural size, i.e. of order one. Using the current knowl- 
edge of the quark mass ratio [40j 



0.47 ±0.08, 



IV. DECAY WIDTHS 

As shown in the last section, the charmed meson loops 
can be neglected. Hence, the LO decay amplitudes are 
given by Eq. (fT0|). Then the decay widths are 



-^rB (md - m u ) 

fir 



—B (m d - m u ) 



(22) 



The pion momenta qi where introduced in the previous 
section. The ratio of these two widths is free of any 
parameter 



= 5.86 ±0.94, (23) 



where the 15% uncertainty comes from neglecting higher 
order terms in either the heavy quark expansion or the 
chiral expansion 



O 



A 



QCD 



0(^1-15%, 



where A x ~ 1 GeV, as well as heavy meson loops. Using 
the experimental value of T (i// — > /i c 7r°) , we predict the 
width of the rj c — > Xco 71 " as 



T (v'c -> Xcoti- ) = 1.5 ± 0.3 ± 0.2 keV, 



(24) 



where the first uncertainty is experimental and the sec- 
ond theoretical due to neglecting higher orders. With 



(27) 



and the LO relation between the pion mass and the quark 
masses m 2 — Bo(m u + md) (neglecting strong isospin 
violation), we get 



B Q (m d ~m u ) 



1 

1 + r 



m ^ = (6.6±2.0) x 10-* GeV, 

(28) 

where the uncertainty is dominated by that of the quark 
mass ratio. Using the value given in Eq. (|28|) . the width 
for the ip' — > h c ir a is 

T (if)' -> h c 7T°) = (0.9 ± 0.6)C 2 keV. (29) 

Since C is of order 1, the above result agrees with the 
BES-III measurement 0.26 ± 0.05 keV well. The agree- 
ment in turn supports the tree-level dominance argued 
for in this paper. It is also consistent with the QCDME 
results reported in Refs. [241 [25|. 



V. DISCUSSION AND SUMMARY 

We demonstrated that, based on a non-relativistic ef- 
fective field theory (NREFT), charmed meson loops are 
highly suppressed by a factor vq 2 / (M^v 2 ) 2 < 1 in 
ip' — > h c n and rj' c —> XcO 71 " ! which are transitions be- 
tween one P-wave and one 5-wave charmonia. The rea- 
son for the suppression is that, due to the small phase 
spaces of the two transitions, the pion momentum is 
much smaller than the approximate kinetic energy of 
the intermediate charmed mesons. The situation is com- 
pletely different for the transitions between two S-wave 
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charmonia. For the ifi' —> J / ifur (tj) , the charmed me- 
son loops are enhanced by a factor of 1 jv compared with 
the tree level contribution (ljjj . There is no factor pro- 
portional to q 2 , and hence the relative size of the pion 
momentum to the kinetic energy of the virtual charmed 
mesons does not have an impact. The difference in these 
two cases is a consequence of the difference in the quan- 
tum numbers of the J/ip and h c which determine their 
coupling to the charmed mesons. 

We note in passing that it follows from our find- 
ings that approximating the ip'h c ?l coupling by that for 
ip'J/ifjr], as done in Ref. 20], is not justified. In the light 
of this one then understands why the estimate made there 
is much larger than the measurement. 

In summary, we have shown that intermediate charmed 
meson loops are highly suppressed in the decays ip 1 — ¥ 
h c ir° and rj' c — > XcO^°, which is completely different from 
the situation of the ip' — > J / iptr (v) ■ We confirmed the 
general power counting arguments given by explicit cal- 
culations both within the NREFT as well as an effec- 
tive Lagrangian approach. By constructing the LO chi- 
ral Lagrangian for the decays, and employing the ex- 
perimental result for T(ip' — ¥ h c ir°), we give a model- 
independent prediction for the width of the r]' c — > Xcqk , 
which can be measured at PANDA, T(rj' c — > XcO n °) — 
1.5 ± 0.4 keV, where the experimental and theoretical 
uncertainties have been summed in quadrature. Note, 
were the transitions dominated by the heavy meson loop 



contributions, the predicted partial decay width would be 
larger by a factor of about 8. An experimental confirma- 
tion of our prediction would provide a strong support for 
the NREFT employed. Having available an effective field 
theory that allows one to study both direct transitions 
as well as those mediated via heavy loops is an impor- 
tant step towards a detailed theoretical understanding of 
charmonium states. 
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